End-to-End Distribution Function of Two-Dimensional Stiff Polymers 

for all Persistence Lengths 
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We set up and solve a recursion relation for all even moments of a two-dimensional stiff polymer 
(Porod-Kratky wormlike chain) and determine from these moments a simple analytic expression for 
the end-to-end distribution at all persistence lengths. 
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I. INTRODUCTION 

In a recent note [l| , two of us set up a recursion rela- 
tion for the even moments of the end-to-end distribution 
of stiff polymers in D dimensions and used the resulting 
moments of high order to construct an simple analytic 
distribution function of the end-to-end distance R. For 
large persistence length £, the result agrees well with per- 
turbative and Monte Carlo results of Wilhelm and Frey 
2] , for small £ with the Daniels correction . 

Recently, Dahr and Chaudhuri |4| have pointed out the 
existence of an interesting dip-structure in the distribu- 
tion function for two dimensions appearing at interme- 
diate ^-values in the spatial density function P^(R) = 
P L {R)/R normalized to unity as / d 2 xP L (R) = 1. The 
often-plotted radial density function Pj,(i?) with the nor- 
malization J Pt, (R) dR — 1 makes this dip almost invisi- 
ble. If one wants this dip tp show up in a simple analytic 
approximation of the type in Ref. 0, several more pa- 
rameters are needed. The purpose of this note is to con- 
struct such an analytic expression which fits excellently 
high-precision Monte Carlo data. 

The end-to-end distribution function of a stiff polymer 
in two dimensions is given by the path integral Q 

P L (H)orJ d(j) b d(j) a V<j>{s)e' Eb/kBT 8 {2) (ft--J o dsu(s)J 

(1.1) 



with the bending energy 



E h 



ds [u'(s)f 



(1.2) 



where u(s) = (cos (j)(s), sin (j)(s)) are the direction vectors 
of the polymer links, and k is the stiffness which defines 
the persistence length £ = 2K/fcsT. Due to the presence 
of the (5^ 2 •'-function in the integrand, the path integral is 
not exactly solvable. It is, however, easy to find arbitrar- 
ily high even moments of the radial distribution of the 
end-to-end distance Pl(R): 



(R 2n ) = J dRR 2n P L (R). 



(1.3) 



This was explained for D-dimensions in Ref. Q. In two 
dimensions, the moments (ll.3|) can be obtained from the 
coefficient of A 2 ™/2 2n (2n)! in the expansion of the integral 



/(r;A)= / d0V(0,T;A), 
Jo 



(1.4) 



in powers of A, evaluated at the euclidean time r = £. 
Here we use natural units with L = 1. The wave func- 
tion ip(6, t; A) is a solution of the Schrodinger equation 
in euclidean time 



d 

d7 



Mr; A), 



where 



H = --A + -Xz 
2 2 



Id 2 1 
2d^ + 2 AC0S ' 



(1.5) 



(1.6) 



II. RECURSIVE SOLUTION OF THE 
SCHROEDINGER EQUATION. 



The function f{L; A) has a spectral representation 
^_^ ^d0^(0)c^(-E^L) ^)(0) 



1=0 



where the ip^> (9) are arbitrarily normalized eigensolu- 
tions of the timc-indcpcndcnt interacting Schrodinger 
equation Htp^'(9) = E^'ipy' {&). Applying perturbation 
theory to this problem, we start from the eigenstates 
|Z) of the unperturbed Hamiltonian Hq = —A/2, with 
eigenvalues E^ = l 2 /2, and form the eigenfunctions 



tp( l >(6) = (6\l) with the explicit form tp(°>{6) = 1/V4tt 
and tp( l '(Q) = cos(l9) / y/n . Note that the ground state 
wave function is not normalized to unity on purpose, for 
later convenience. Now we set up a recursion scheme for 
the expansion coefficients 7^ and ep of the eigenfunc- 
tions and their energies: 



,(0\ - 



= £ 7 gA* |0, ef A'. (2.2) 



l',i=0 



3=0 



2 



The procedure is described in 0, |(| . The properties of 
the unperturbed system for A = determine the initial 
conditions for the recursion: 



(0 x 



(9 X 



So) 



f/2. (2.3) 



To proceed, we need the the matrix elements of the per- 
turbing Hamiltonian Hi in the basis of the unperturbed 
eigenstates. They are simply (n\Hj\n±l) = A/2. Insert- 
ing the expansions (|2.2|l into the Schrodinger equation 
1)1. 5[1 . projecting the result onto some base vector (k\, 
and extracting the coefhcient of A 1 , we obtain the follow- 
ing recursion relations: 



and 



(o) 



7i,t-D 



(0 



r- 



(V) 



(7 i ( | ) M _ 1 +7^ M _ 1 )/2, (2.4) 



i=i 



(0 (0 
: j 7o,<- 



(2.5) 



7fc 



(9 



1 



P-k 2 



'fe+l,t— 1 T 'fc-l,*-l A^i 'k,i-j 



(2.6) 

Starting from the initial values (|2.3|) , these recursion rela- 
tions determine successively the higher-order expansion 
coefficients in Q2.2p . Inserting the resulting expansions 
(|2.2|) into Eq. I|2.1|l . only the constant parts in <p^(9) 
which are independent of 9 will survive the integration 
in the numerators. Therefore tp® (9) in the numerators 
of 12.1fl may be replaced by the constants: 



rsymm 



r 27T OO 



A' 



(2.7) 



the factor 1/2 reflecting the special normalization of 



,(0) 



(9). The denominators of (|2.1|) become explicitly 



Jd9^(9) ^H^) = E (l^| 2 /2 + El^il 2 ) A 



1 - 1 y2i 

(2.8) 



where the sum over i is limited by the power of A 2 up 
to which we want to carry the perturbation series; also 
I' is restricted to a finite number of terms, because of a 
band-diagonal structure of the jj-P, (see Q). Extracting 
the coefficients of the power expansion in A from (|2.1() 
we obtain all desired moments of the end-to-end distri- 
bution, the lowest two being: 



(R 2 



(R 4 



'fa 



i 



:8L 2 £ 2 -L£ 3 [ :-!() 
'87 392 



-L/i 



40 



(2.9) 



-ALU 



9 



18 



(2.10) 



The calculation of higher moments is straightforward 
with a Mathematica program, which we have placed on 
the internet in notebook form ML The above lowest mo- 
ments agree with those in Ref. 



III. END-TO-END DISTRIBUTION 

As in the previous paper, we shall now set up an an- 
alytic distribution function for Pi(r) = P^(R), where 
r '= R/L: 



Pl{t) = (a + a 2 £ +a 4 :c +a 6 cc b )ar(l - x p ) 



/3\m 



(3.1) 



The parameters ao, . . . , a^, fc, (3, m are functions of £/L 
and arc determined by forcing the moments of (|3.1|l . 



2n\ 



- 2n 



r 2n+1 p L (r) dr, 



(3.2) 



to fit the exact moments in the range < n < 
Max(6, 1Q£/L). For floppy chains with £ < 1 we set 
k = 0; for very floppy chains with £ < 1/20 the choice 
a2 = a4 = a6 = guarantees excellent accuracy. A 
comparison of Pl(t) with Monte-Carlo data is shown in 
Fig. ^ an d with higher resolution near R = in Fig. 
The associated coefficients arc listed in Table |l| The 



10 




0.2 



0.4 



0.8 



0.6 

r = R/L 

FIG. 1: End-to-end distribution Pl(t) in D — 2 dimensions 
as a function of r = R/L for various values of the stiffness 
i/L = 1/15, 1/10, 1/5, 1/4, 3/10, 7/20, 2/5, 1/2, 1, 2. The 
solid curves show the model functions 13. II with parameters 
from Table U The dots show the Monte Carlo data. 



calculation of the coefficients in l|3.1|l requires some care 
to guarantee sensitivity to possible local minima, and to 
avoid running into unphysical oscillations. The latter 
may arise from the existence of polynomials in which all 
moments lower than some n vanish. Such oscillations are 
avoided by controlling the high moments and using only 
low polynomial coefficients in (|3.1|l . For £ > 1/20, a more 
involved strategy is necessary to avoid low-quality local 
solutions. We proceed as follows: 



3 



£ 


do 


ai 


CE4 


Ct6 


k 


m 


P 


1/400 


400.0 














196.784 


1.99496 


1/100 


100.0 














47.5378 


1.98197 


1/50 


50.0 














22.9930 


1.97564 


1/30 


29.5302 


-58.9195 


77.9373 


-87.3526 





12.0224 


2.00505 


1/15 


14.0952 


-29.8504 


66.8842 


-68.1985 





5.62896 


2.21169 


1/10 


9.20629 


-34.7515 


50.6223 


-26.2289 





13.2737 


10.5486 


1/5 


4.18239 


-7.45808 


11.616 


-7.30855 





10.2031 


16.6444 


1/4 


3.12655 


-4.9930 


13.1086 


-10.0222 





9.42195 


20.0750 


3/10 


2.38054 


-3.38168 


12.8823 


-9.51483 





9.16782 


23.0164 


7/20 


1.82132 


-2.062292 


11.2343 


-7.24306 





8.84230 


25.5206 


2/5 


1.39171 


-0.952158 


8.94986 


-4.33545 





8.49552 


27.9120 


1/2 


0.800939 


0.647524 


4.36711 


1.36933 





8.06681 


33.2814 


1 


42.8376 


-173.308 


263.327 


-123.515 


4.9880 


11.4933 


85.8428 


2 


504.624 


-1832.52 


2271.51 


-925.829 


13.4792 


30.4949 


244.143 



TABLE I: Coefficients of the analytic distribution function for in Eq. 13. H for various values of the persistence length £. 

They are obtained by making six or seven even moments of Pl(t) agree with the exact ones. 





log p L (r) 


i/L = l/40o\ 
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FIG. 2: End-to-end distribution near R = in D = 2 di- 
mensions plotted logarithmically as a function of r — R/L 
for various values of the stiffness for floppy polymers with 
£/L = 1/400, 1/100, 1/50, 1/30. The solid curves show the 
model functions (13. II with parameters from Table Q The dots 
represent Monte Carlo data. 



• In a first step we set 02 = 04 = a& = 0, (3 = 2, 
and determine preliminary values for k and m by 
fitting two higher moments of n near 10 £/Z. The 
first coefficient ao is fixed by normalization. This 
gives a reasonable starting value for m. 

• In a second step, we introduce one more of the 
higher moments and improve the solution for k, m, 
and (3. 

• Next we solve for the coefficients aj by bringing 
yet more moments into play. If £/Z < 1, we take 



k = and solve for the Oj, keeping /3 and m fixed, 
based on four properly chosen moments. Then we 
solve for j3 and m keeping the a,j fixed, based on a 
choice of two moments. This alternating procedure 
is repeated three times. Finally, we solve for the 
cij , f3, and m simultaneously, based on six properly 
chosen moments. 

For £/L > 1, we proceed similarly, but allow for k 7^ 
0. The search for good coefficients aj alternating 
with a search for good k, f3, and m is repeated until 
it converges, with no further attempt to solve once 
more for all seven parameters simultaneously. 



We check the quality of our simple distribution func- 
tion l|3.1fl with the parameters of Table [I] by calculating 
its moments and comparing them with the exact ones. 
The comparison is shown in Table ITT1 for a large range of 
the persistence length As a measure of the quality of 
approximation we use the quantity S, listed in the second 
column of Table [n] which sums up all squared deviations 
of the moments of the model from the exact ones in a 
relevant range of £: 



£({) 



\ 



JV 

71 = 



(3.3) 



where we have extended the sum up to the moments of 
order N = 12 for f/Z < 1/5, and N = 24 for f/Z > 1/5. 

Let us also convince ourselves of the high accuracy of 
our Monte Carlo data for the end-to-end distribution in 
Figs. Hand [21 by listing the maximal deviation 



^abs 



sup 

71=0 



\(R 



2h \ 



IMC 



(R 



2n\ 



(3.4) 



4 





E 


A a bs 




1/400 


3 x 10~ 12 


0.000 126 


0.9%(16) 


1/100 


2 x 10" 13 


0.000 033 


8%(16) 


1/50 


1 x 10~ 10 


0.000 071 


1%(16) 


1/30 


5 x 10~ 9 


0.000 717 


8%(16) 


1/15 


2 x 10~ e 


0.000 057 


0.9%(16) 


1/10 


5 x 10~ 5 


0.000 038 


0.8%(24) 


1/5 


4 x 10~ 5 


0.000 048 


0.5% (24) 


1/4 


9 x 10 -5 


0.000 048 


0.4%(24) 


3/10 


13 x 10~ 5 


0.000 047 


0.3%(24) 


7/20 


2 x 10~ 4 


0.000 087 


0.8%(36) 


2/5 


2 x 10~ 4 


0.000 100 


0.8%(36) 


1/2 


2 x 10~ 4 


0.000 139 


0.8%(36) 


1 


2 x 10~ 4 


0.000 217 


0.4%(48) 


2 


8 x 10" 5 


0.004 705 


1.6%(48) 



TABLE II: To illustrate the accuracy of our analytic approx- 
imation 13.11 of the end-to-end distribution we list the quan- 
tity E of Eq. £S) which measures the deviation of the even 
moments from the exact ones. The other two columns il- 
lustrate the accuracy of the Monte Carlo data for the end- 
to-end distribution by listing the maximal deviation A a b s of 
its moments and the relative deviation A re i(iV max ) up to the 
moment iV m a X . 



of its moments and the relative deviation 

A re i(AU x ) = sup | (R 2n ) MC /(R 2n ) PL ) - 1| (3.5) 

n=0 



up to the moment Aj 



Remarkably, in spite of the simplicity of the model, it is 
a nontrivial task to obtain accurate Monte Carlo results 
for Pj,(R) near R = which are sensitively displayed in 
the plots of Fig. but which are almost ignored by the 
moments. The reason for the difficulty is the small con- 
figuration space for the small-R data since the binning 
of the data is done in R to estimate the density Pl(R). 
One is caught in the competition between large system- 
atic errors resulting from a too large bin size AR, and 
statistical errors from a too small AR. As a compro- 
mise we employed in our simulations a uniform bin size 
AR/ L = 0.01 which in combination with a single-cluster 
update procedure and a statistics of 10 8 sampled chains 
yielded satisfactory accuracy near R = 0, as shown in 
the high-resolution plot in Fig. 
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